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Abstract. We review the current status of recombination and coalescence 
models that have been successfully applied to describe hadronization in heavy 
ion collisions at RHIC energies. Basic concepts as well as actual implementations 
of the idea are discussed. We try to evaluate where we stand in our understanding 
at the moment and what remains to be done in the future. 
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1. Introduction 

Since the last Quark Matter conference, held in 2002, the idea of quark recombination 
or coalescence (ReCo) had a very rapid and successful career as a model for 
hadronization. Recombination of quarks was first formulated over two decades ago 
by R. C. Hwa and collaborators as a model for hadronization in the fragmentation 
region of hadron-hadron collisions . There were early efforts to connect it to other 
descriptions of hadronization like fragmentation [21 El- For heavy ion collisions, quark 
coalescence ideas first appeared soon after that 0] and were later on successfully used 
in ALCOR 0. 

The recent revival of recombination/coalescence began when it was realized [H] 
that the elliptic flow z; 2 measured at RHIC obeys a simple valence quark scaling, that 
naturally arises from a recombination picture. Soon after that, it was pointed out 
[71 [8] that other RHIC puzzles, like the anomalous enhancement of baryons and the 
absence of nuclear suppression in baryon spectra can be explained by ReCo as well. 

This talk seeks to review the current status of recombination/coalescence models. 
We will revisit the fundamental concepts and experimental findings that support 
recombination. We will then discuss different implementations of the model and its 
limitations. We close with an outlook on future developments. 



2. Basic Ideas 



Hadronization has always been a very difficult aspect of strong interaction processes. 
A lot of effort went into the invention of methods to work around this non-perturbative 
phenomenon, with the effect that our knowledge of hadronization dynamics is still very 
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limited. One example is the formulation of fragmentation or "quark decay" functions 
which rely on the common concept of separating long and short distance dynamics 
UJ- The cross section of hadron production in e + e~, lepton-hadron or hadron-hadron 
collisions can be written as a convolution dan = do~ a C§) D a —,u of a parton production 
cross section da a with a fragmentation function D a ^n(z, n) ^Hl- D a ^n{z, n) gives 
the probability to find a hadron H in the hadronizing parton a with a momentum 
fraction z, < z < 1. /i> Aqcd is a perturbative scale that is set by the transverse 
momentum. 

The z-dependencc of fragmentation functions is not calculable in perturbation 
theory. However, they are universal, i.e. process independent, objects. Fragmentation 
functions for the most important hadrons have been measured, mainly in e + e~ 
collisions, and are available in parametrized form. Nevertheless, the availability of 
measured fragmentation functions should not lead to the impression that we have a 
good understanding of the underlying hadronization process. 

One important question is the following: at which scales (i.e. in which Pt region) 
is leading twist perturbative QCD (pQCD) reliable? It has been shown that next- 
to-leading order (NLO) calculations of ir° production in pp collisions at RHIC, using 
modern parametrizations of fragmentation functions are in good agreement with 
PHENIX data down to surprisingly low pion Pt of 1 GeV/c |12|. 

On the other hand, there are clear signs that the fragmentation concept is not 
working at very low Pt- A good example is the leading particle effect. In the very 
forward (and low Pt) region of hadron-hadron collisions pQ, the composition of particle 
species deviates from expectations in a fragmentation picture. This is impressively 
confirmed by recent experimental results, e.g. from the FNAL E791 collaboration 
|13j . With a 7r~ beam impinging on a fixed nuclear target, they measure a D~ /D + 
asymmetry that goes to 1 in the very forward direction. While fragmentation would 
predict this asymmetry to be very close to 0, one can understand this effect starting 
from a recombining cc pair produced in the collision. The recombination cd —> D~ is 
enhanced with respect to cd — > D + , because the d is a valence quark in the beam n~ 
while the d is only a sea quark ^1] . 

Other examples are provided by the RHIC experiments. A proton/pion ratio ~ 1 
was observed in central Au+Au collisions between Pt = 1.5 and 4 GeV/c [T^J. This is 
in contradiction to pQCD calculations, that give p/ir ~ 0.1 . . . 0.2 jJSj. Similar results 
for A/K® ^2] suggest that there is a general pattern of baryon enhancement at RHIC 
energies. This enhancement is so strong that it neutralizes the strong jet quenching 
observed for mesons at RHIC |18|. The nuclear modification factor Raa for baryons 
is close to 1 up to 4 GeV/c ^3 EI EH- Thus, although pQCD seems to work very 
well in p +p for P T above 1 GeV/c, this is not the case in Au+Au even at 4 GeV/c. 

The two examples have in common that hadronization takes place in a phase 
space filled with partons, either beam remnants or the hot and dense medium created 
in heavy ion collisions. This is very different from e + e~ collisions where phase space 
is nearly empty. One can include corrections to single parton fragmentation in terms 
of higher twist fragmentation. Such contributions include the process of two or more 
partons fragmenting into hadrons. However, nothing is known about these higher twist 
corrections. Instead of the rather complicated twist expansion, let us directly look at 
the limiting case of a phase space densely packed with partons. What will happen 
upon hadronization? The most simple picture is that the quarks and antiquarks 
that constitute the valence quark structure of a hadron (having the correct quantum 
numbers), recombine/coalesce together to form this hadron. 
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Figure 1. Recombination and fragmentation for a meson at Pt ~ 6 GeV/c, 
starting from a parton spectrum with steep slope (solid line). Fragmentation 
requires a single parton with transverse momentum larger than to start with, 
while recombination is possible with two partons that have roughly Pt/2 each. 
The competition between both processes is decided by slope and normalization of 
the parton spectrum. 



In fragmentation the entire parton content of the hadrons has to come from 
gluons and qq pairs emitted from the fragmenting parton inside the "black box" called 
fragmentation function. For recombination, it is assumed that no further branching 
of partons occurs. Apparently these are limiting cases. There must be a smooth 
transition between both extremes as a function of phase space density. On the other 
hand, one might speculate that the absence of additional branching is a hint that some 
sort of equilibrium is reached in the parton phase. Figure^shows how both processes 
form a meson. 

Which partons enter the recombination process? In most implementations, quarks 
are assumed to be non-perturbative and to have an effective mass. This, and the fact 
that only valence quarks are involved, leads to the interpretation of these degrees 
of freedom as constituent quarks. Gluons do not participate in recombination and 
are used to dress the quarks. The alternative to the constituent quark picture is to 
explicitly convert gluons into qq pairs. 

3. Mathematical Formulation 

ReCo can be formulated in terms of Wigner functions. The yield of mesons M 
coalescing from two partons a, b is given by 

dN M ^ f d 3 R d 3 qd 3 r^ r ( ' r P n r P \ . . 

a.b 

W a b and $m are the Wigner functions of the partons and the meson respectively, P 
and R are the momentum and spatial coordinate of the meson and the sum runs over 
all possible combinations of quantum numbers, essentially leading to a degeneracy 
factor Cm- The generalization of this formula for baryons is straightforward [TJ]. The 
Wigner function for the partons is usually factorized into classical one-particle phase 
space distributions, W a b(faiPai r biPb) = w a {r a) p a )wb(rb,pi,)- This assumes that the 
partons are completely uncorrelated before hadronization. We will come back to this 
point later. 
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Figure 2. Elliptic flow V2 scaled with the number n of valence quarks vs Pr/n 
for four different hadrons. Data are taken from STAR (K°, A + A) [l^l and 
PHENIX (7r+,p) El- 



We can immediately draw some conclusions. Suppose the parton spectrum is 
exponential in transverse momentum px, w = Ae~ PT ^ T ', with some slope parameter 
1/T. Then fragmentation and recombination would provide meson spectra 

dN llag /d 2 P T D oc Ae- PT/{Z)T (D), dN rcco /d 2 P T oc A 2 e" PT/T (2) 

respectively. Here (D) and (z) are average values of the fragmentation function and the 
scaling variable z. Since (z) < 1, fragmentation is less effective than recombination 
on an exponential spectrum, as long as the normalization A is not too small. On 
the other hand, if the parton spectrum has power law form, w = Bp^ a , the yields 
are dNf rajg /d 2 PT ~ P^ a and dN Teco /d 2 PT ~ P^ 2a for mesons. This implies that 
fragmentation will dominate at large Pp for power law spectra (which are predicted 
by perturbative QCD). 

If the parton spectrum is exponential, one can check that with increasing Pt the 
result is less and less sensitive to the momentum dependence of the hadronic Wigner 
function (i.e. the shape of the momentum space wave functions used to model it). In 
fact, it seems that for Pt > 2 GeV/c one can safely take the momentum space wave 
function to have zero width (i.e. to be a 6- function) |16| . 

Let us assume that the parton phase exhibits an azimuthal anisotropy. The elliptic 
component of this asymmetry is described by the elliptic flow coefficient V2(pt) |2(Jj. 
ReCo predicts the resulting elliptic flow for hadrons to be El 123 

v^(P T ) = 2v 2 (^j , v ^p T ) = 3v 2 (^j (3) 

for mesons and baryons respectively JSj. One should note that the derivation of 
these scaling laws uses narrow, (5-shaped wave functions. See ^Sj for the case of wave 
functions of finite width. As can seen be in Figure the scaling laws are impressively 
confirmed by experimental data for Pt > 1 GeV/c, while at lower Pt the mass 
of the hadron determines the elliptic flow, well described by hydrodynamics |191 122] . 
Deviations for pions are discussed below. Above ~ 5 GeV/c a perturbative mechanism, 
driven by jet quenching, should take over from recombination, but experimental data 
are not yet conclusive |lfi| . 
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4. Implementations 

The development of ReCo in the past 18 months was mainly driven by four groups who 
published work using their individual implementations of the ReCo concept. These 
groups are (in parentheses the short form used thereafter) Duke/Minnesota/Kyoto 
(Duke) Ohio State/Wayne State (Ohio) 1121211157], Oregon (2H1ES1 

and Texas A& M/Budapest (TAM) |H EDI EH E2 

The Duke and Oregon groups try to evaluate Equation Q analytically using 
certain assumptions that essentially boil it down to a convolution of one-particle 
distributions and hadron wave functions in longitudinal momentum space (longitudinal 
with respect to the hadron momentum). As an example, the Duke group writes the 
meson spectrum from recombination as |16| 



where da integrates over the hadronization hypersurface E, u is the four vector 
orthogonal to E, x is the momentum fraction of parton a in the meson, P + is the 
light cone momentum and §m the wave function of the meson. The parton phase 
that undergoes hadronization is assumed to have an exponential part at low px (soft 
partons) and a power law tail at high px (hard partons). For central Au+Au collisions 
at RHIC, the Duke group uses a thermal distribution with temperature T = 175 MeV 
and average radial flow velocity v = 0.55c for the soft partons, while the hard partons 
are taken from a minijet calculation jSH] including energy loss |16| . 

The Ohio and TAM groups developed Monte Carlo implementations of the 
recombination process. These can be connected to string or parton cascade models 
that prepare the parton state before hadronization. One main difference lies in the 
treatment of the connection between soft and hard partons. While the Duke group 
strictly separates soft and hard physics, allowing only the soft partons to recombinc 
and only the hard partons to fragment, the TAM group includes additional coalescence 
of soft and hard partons [5] |3U|. The Oregon group carries this a step further 
and replaces fragmentation functions by a scenario where minijet partons develop 
a shower which subsequently recombines. In this model, they are able to describe 
fragmentation functions reasonably well. Applied to heavy ion collisions, this allows 
them to introduce recombination of soft partons with shower partons 2!) . 

All four groups describe hadron data from RHIC at intermediate and large 
transverse momenta very well. The calculations for spectra reproduce the noticeable 
kink in the data around Pt = 4 GeV/c for mesons and 6 GeV/c for baryons 
coming from the transition from soft (recombination) to hard (perturbative) particle 
production. Soft-hard coalescence can improve the fit to data points in the transition 
region, which is then extended to even higher Py. Figure [21 shows a result obtained 
by the Duke group compared with RHIC data. The baryon/meson ratios in ReCo are 
essentially given by the ratio of degeneracy factors Cb/Cm, leading naturally to an 
enhanced proton/pion ratio. The most recent results can be found in [T?)l l2~§| . 

As already mentioned, the elliptic flow of all measured hadron species, tt, K, p, A, 
2 and f2, follows the scaling law JIJ for Pt > 1 GeV/c. This permits to unambiguously 
extract V2(j>t) for partons. Surprisingly, the elliptic flow of strange quarks is the same 
as for light quarks Slight deviations from the scaling can be seen for pions. 

This can be traced back to the small mass of the pion and the fact that most pions 
do not hadronize directly but are from p decays |32| . 
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Figure 3. Spectra of 7r°, p, K® and A + A as a function of Pt at midrapidity 
for central Au+Au collisions at ^/s = 200 GeV/c 1161 . RHIC data are taken from 
PHENIX (7r°,p) E3[TS] and STAR (K°,A + A) [T7|. 



5. ReCo Challenges and Answers 

Despite its success, several aspects of ReCo are problematic. For one, recombination 
in its current implementation cannot describe the bulk production of hadronic matter. 
This seems to contradict the claim that ReCo is the right choice for very dense parton 
systems. However, one should note that the ReCo formula is only for the average 
meson created by recombination (like the fragmentation formula is for the average 
hadron produced from one parton). None of these formulas describe the exclusive 
hadron content of the system. Note that in (|TJ the mesons scale with the square of 
the parton density, whereas the total number of mesons has to scale linearly, of course. 
While ReCo as a concept is certainly correct also for the bulk of hadron production, 
the simple formula |Q} does not describe this. 

The second issue is that Equation |JTJ conserves momentum, but not energy. In 
general, it is not possible to have energy conservation in 2 — > 1 and 3 — ► 1 processes. 
In reality, particles in this strongly interacting environment will be off-shell, making 
energy conservation possible. However, instead of taking displacements from the 
mass shell into account, it is more convenient to have particles on the mass shell 
and permit small violations of energy conservation. This is a safe procedure as long as 
the violations are small compared to the total energy E of the hadron, i.e. for above 
~ 2 GeV/c A special role is played by Goldstone bosons. Their description is 
particularly difficult in a picture where the coalescing quarks have constituent masses. 

The situation can be improved by taking into account coalescence of resonances 
and their subsequent decay. This introduces 2^2 and 3^2 processes and energy 
can be exactly conserved. Including the p resonance, decaying into two pions, helps to 
cure the problem that pions are too heavy in a constituent quark picture | 32| . Another 
critical question concerns entropy. Apparently, recombination reduces the degrees of 
freedom in the system, therefore leading to a decrease in entropy. We note that this 
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statement is not relevant as long as we do not address total particle numbers. But in 
any case, including resonances considerably improves this situation as well. 

6. For the Future 

One can think of a long list of problems that should be addressed in the near future. It 
was pointed out that the study of resonances within the recombination model provides 
new insights about interactions in the hadronic phase and could even be used to pin 
down properties of resonances and exotic states [23] . It also has to be decided whether 
charmed hadrons follow the recombination systematics. The elliptic flow of charm 
quarks has to be measured |261 l3"T] . First results on higher harmonics are available 
now |35) . They provide novel tests for the ReCo concept and can be used to further pin 
down the partonic phase Furthermore, the role of soft-hard coalescence has to be 
investigated in more detail, in particular with respect to hadron-hadron correlations. 

RHIC data indicates that strong hadron-hadron correlations are present at 
intermediate Pt, where ReCo dominates [37]. Such correlations are expected for 
the fragmentation process where several hadrons emerge from the same parton. It 
was argued that hadrons from recombination are emitted statistically so that no 
correlations should be observable. This is not true. Correlations in the parton phase 
are naturally translated into correlations in the hadron phase by the recombination 
process. We immediately conclude that there must be non-negligible correlations 
between soft partons. 

Let us assume we want to describe the recombination of four partons 1,2,3,4 into 
two mesons A, B. In analogy to Equation QJ, a 4-parton Wigner function VF1234 
is needed to describe double meson production. We can include correlations in the 
parton phase by replacing the simple single particle factorization by 



where CV, is a 2-parton correlation function and N is a normalization factor. This 
will lead to non trivial correlation functions between mesons. The correlations in the 
parton phase could originate from interactions between hard and soft partons and 
would be naturally present in soft-hard recombination. It remains to be seen whether 
the correlations measured at RHIC will consistently fit into this picture. But if this 
is the case, ReCo will provide a fascinating new picture of the partonic phase. 

What will happen at LHC? Part of the ReCo success story is, that jet quenching 
is so strong at RHIC, suppressing hard processes in the final state by a factor of 4. 
Hard processes will be more abundant at LHC, but increased jet quenching and a 
brighter thermal source will probably overwhelm them up to even higher Pt than at 



7. Summary 

Quark recombination/coalescence is a successful model to describe hadronization in 
dense parton systems. Central Au+Au collisions provide a partonic medium that 
is sufficiently dense for coalescence of soft partons to overcome fragmentation at 
intermediate transverse momenta between 2 and 5 GeV/c. Soft-hard coalescence could 
push this region to even higher Pt- 




(5) 
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